OcBSA: An NGS-based bulk segregant analysis tool for outcross populations.

Constructing inbred lines for self-incompatible species and species with long generation times is challenging, making the use of F1 outcross/segregating populations the main strategy for genetic studies of such species. However, there is a lack of dedicated algorithms/tools for rapid quantitative trait locus (QTL) mapping using the F1 populations. To this end, we have designed and developed an algorithm/tool called OcBSA specifically for QTL mapping of F1 populations. OcBSA transforms the four-haplotype inheritance problem from the two heterozygous diploid parents of the F1 population into the two-haplotype inheritance problem common in current genetic studies by removing the two haplotypes from the heterozygous parent that do not contribute to phenotype segregation in the F1 population. Testing of OcBSA on 1800 simulated F1 populations demonstrated its advantages over other currently available tools in terms of sensitivity and accuracy. In addition, the broad applicability of OcBSA was validated by QTL mapping using seven reported F1 populations of apple, pear, peach, citrus, grape, tea, and rice. We also used OcBSA to map the QTL for flower color in a newly constructed F1 population of potato generated in this study. The OcBSA mapping result was verified by the insertion or deletion markers to be consistent with a previously reported locus harboring the ANTHOCYANIN 2 gene, which regulates potato flower color. Taken together, these results highlight the power and broad utility of OcBSA for QTL mapping using F1 populations and thus a great potential for functional gene mining in outcrossing species. For ease of use, we have developed both Windows and Linux versions of OcBSA, which are freely available at: https://gitee.com/Bioinformaticslab/OcBSA.

Characterization of pollen tube development in distant hybridization of Chinese cork oak (Quercus variabilis L.).

MAIN CONCLUSION: This work mainly found that the stigma and style of Q. variabilis did not completely lose the specific recognition towards heterologous pollen, a fact which is different from previous studies. Quercus is the foundation species in the Northern Hemisphere, with extreme prevalence for interspecific hybridization. It is not yet entirely understood whether or how the pollen tube-female tissue interaction contributes to the "extensive hybridization" in oaks. Pollen storage conditions correlate with distant hybridization. We conducted hybridization experiments with Q. variabilis as female and Q. variabilis and Q. mongolica as male parents. And the differences in pollen tube (PT) development between intra- and distant interspecific hybridization were studied by fluorescence microscopy and scanning electron microscopy (SEM). Our results showed that -20 °C allowed pollen of both species to maintain some viability. Both Q. variabilis and Q. mongolica pollen germinated profusely on the stigmas. SEM results indicated that in the intraspecific hybridization, Q. variabilis pollen started to germinate at 6 h after pollination (hap), PTs elongated significantly at 12 hap, and entered the stigma at 24 hap. By contrast, Q. mongolica pollen germinated at 15 hap, and the PTs entered the stigma at 27 hap. By fluorescence microscopical studies it was observed that some PTs of Q. variabilis gathered at the style-joining at 96 hap, unlike the Q. mongolica which reached the style junction at 144 hap. The above results indicate that the abundant germination of heterologous pollen (HP) on the stigma and the "Feeble specificity recognition" of the stigma and transmitting tract to HP may create opportunities for the "extensive hybridization" of oaks. This work provides a sexual developmental reference for clarifying the causes of Quercus "extensive hybridization".

Intermittent fasting and caloric restriction interact with genetics to shape physiological health in mice.

Dietary interventions can dramatically affect physiological health and organismal lifespan. The degree to which organismal health is improved depends upon genotype and the severity of dietary intervention, but neither the effects of these factors, nor their interaction, have been quantified in an outbred population. Moreover, it is not well understood what physiological changes occur shortly after dietary change and how these may affect the health of an adult population. In this article, we investigated the effect of 6-month exposure of either caloric restriction (CR) or intermittent fasting (IF) on a broad range of physiological traits in 960 1-year old Diversity Outbred mice. We found CR and IF affected distinct aspects of physiology and neither the magnitude nor the direction (beneficial or detrimental) of effects were concordant with the severity of the intervention. In addition to the effects of diet, genetic variation significantly affected 31 of 36 traits (heritabilities ranged from 0.04 to 0.65). We observed significant covariation between many traits that was due to both diet and genetics and quantified these effects with phenotypic and genetic correlations. We genetically mapped 16 diet-independent and 2 diet-dependent significant quantitative trait loci, both of which were associated with cardiac physiology. Collectively, these results demonstrate the degree to which diet and genetics interact to shape the physiological health of adult mice following 6 months of dietary intervention.

Relating multivariate shapes to genescapes using phenotype-biological process associations for craniofacial shape.

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high-dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait - craniofacial shape - in a sample of n = 1145 mice from the Diversity Outbred (DO) experimental line. We use biological process Gene Ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (process MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application that provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex trait genetics.

Genome-wide association study using diversity outcross mice identified candidate genes of pancreatic cancer.

To understand the genetic causes of pancreatic cancer (PC), we conducted a genome-wide association study (GWAS) using the diversity outbred (DO) mice population to identify susceptibility genes underlying 7,12-dimethylbenzanthraene (DMBA) induced PC. The phenotype studied was the percent PC lesion area in the DO mice population. We genotyped 7851 SNP markers specifically designed for DO mice across the whole mouse genome. Four susceptibility genes with P values exceeding the genome-wide threshold for percent PC lesion area (P < 2.37 × 10-6) were identified, i.e., Epha4, Gpc5, Kcnj6, Arid1b. The most significant SNP of Gpc5 (UNC140360310) that is associated with PC lesion area in mice also significantly influences the Gpc5 expression, suggesting that this Gpc5 SNP exerts its role in PC through cis-regulating the gene expression of Gpc5. Together, our data supported that Gpc5 as a tumor suppressor gene involved in the etiology of PC.

Development and characterization of transgenic dominant male sterile rice toward an outcross-based breeding system.

Genomic selection is attracting attention in the field of crop breeding. To apply genomic selection effectively for autogamous (self-pollinating) crops, an efficient outcross system is desired. Since dominant male sterility is a powerful tool for easy and successive outcross of autogamous crops, we developed transgenic dominant male sterile rice (Oryza sativa L.) using the barnase gene that is expressed by the tapetum-specific promoter BoA9. Barnase-induced male sterile rice No. 10 (BMS10) was selected for its stable male sterility and normal growth characteristics. The BMS10 flowering habits, including heading date, flowering date, and daily flowering time of BMS10 tended to be delayed compared to wild type. When BMS10 and wild type were placed side-by-side and crossed under an open-pollinating condition, the seed-setting rate was <1.5%. When the clipping method was used to avoid the influence of late flowering habits, the seed-setting rate of BMS10 increased to a maximum of 86.4%. Although flowering synchronicity should be improved to increase the seed-setting rate, our results showed that this system can produce stable transgenic male sterility with normal female fertility in rice. The transgenic male sterile rice would promote a genomic selection-based breeding system in rice.

Resources for Systems Genetics.

A key characteristic of systems genetics is its reliance on populations that vary to a greater or lesser degree in genetic complexity-from highly admixed populations such as the Collaborative Cross and Diversity Outcross to relatively simple crosses such as sets of consomic strains and reduced complexity crosses. This protocol is intended to help investigators make more informed decisions about choices of resources given different types of questions. We consider factors such as costs, availability, and ease of breeding for common scenarios. In general, we recommend using complementary resources and minimizing depth of resampling of any given genome or strain.

Reconstruction of Genome Ancestry Blocks in Multiparental Populations.

We present a general hidden Markov model framework called R: econstructing A: ncestry B: locks BIT: by bit (RABBIT) for reconstructing genome ancestry blocks from single-nucleotide polymorphism (SNP) array data, a required step for quantitative trait locus (QTL) mapping. The framework can be applied to a wide range of mapping populations such as the Arabidopsis multiparent advanced generation intercross (MAGIC), the mouse Collaborative Cross (CC), and the diversity outcross (DO) for both autosomes and X chromosomes if they exist. The model underlying RABBIT accounts for the joint pattern of recombination breakpoints between two homologous chromosomes and missing data and allelic typing errors in the genotype data of both sampled individuals and founders. Studies on simulated data of the MAGIC and the CC and real data of the MAGIC, the DO, and the CC demonstrate that RABBIT is more robust and accurate in reconstructing recombination bin maps than some commonly used methods.

Modeling X-linked ancestral origins in multiparental populations.

The models for the mosaic structure of an individual's genome from multiparental populations have been developed primarily for autosomes, whereas X chromosomes receive very little attention. In this paper, we extend our previous approach to model ancestral origin processes along two X chromosomes in a mapping population, which is necessary for developing hidden Markov models in the reconstruction of ancestry blocks for X-linked quantitative trait locus mapping. The model accounts for the joint recombination pattern, the asymmetry between maternally and paternally derived X chromosomes, and the finiteness of population size. The model can be applied to various mapping populations such as the advanced intercross lines (AIL), the Collaborative Cross (CC), the heterogeneous stock (HS), the Diversity Outcross (DO), and the Drosophila synthetic population resource (DSPR). We further derive the map expansion, density (per Morgan) of recombination breakpoints, in advanced intercross populations with L inbred founders under the limit of an infinitely large population size. The analytic results show that for X chromosomes the genetic map expands linearly at a rate (per generation) of two-thirds times 1 - 10/(9L) for the AIL, and at a rate of two-thirds times 1 - 1/L for the DO and the HS, whereas for autosomes the map expands at a rate of 1 - 1/L for the AIL, the DO, and the HS.

Floral display size, conspecific density and florivory affect fruit set in natural populations of Phlox hirsuta, an endangered species.

BACKGROUND AND AIMS: Natural variation in fruit and seed set may be explained by factors that affect the composition of pollen grains on stigmas. Self-incompatible species require compatible outcross pollen grains to produce seeds. The siring success of outcross pollen grains, however, can be hindered if self (or other incompatible) pollen grains co-occur on stigmas. This study identifies factors that determine fruit set in Phlox hirsuta, a self-sterile endangered species that is prone to self-pollination, and its associated fitness costs. METHODS: Multiple linear regressions were used to identify factors that explain variation in percentage fruit set within three of the five known populations of this endangered species. Florivorous beetle density, petal colour, floral display size, local conspecific density and pre-dispersal seed predation were quantified and their effects on the ability of flowers to produce fruits were assessed. KEY RESULTS: In all three populations, percentage fruit set decreased as florivorous beetle density increased and as floral display size increased. The effect of floral display size on fruit set, however, often depended on the density of nearby conspecific plants. High local conspecific densities offset - even reversed - the negative effects of floral display size on percentage fruit set. Seed predation by mammals decreased fruit set in one population. CONCLUSIONS: The results indicate that seed production in P. hirsuta can be maximized by selectively augmenting populations in areas containing isolated large plants, by reducing the population sizes of florivorous beetles and by excluding mammals that consume unripe fruits.

Saccharomyces cerevisiae: a sexy yeast with a prion problem.

Yeast prions are infectious proteins that spread exclusively by mating. The frequency of prions in the wild therefore largely reflects the rate of spread by mating counterbalanced by prion growth slowing effects in the host. We recently showed that the frequency of outcross mating is about 1% of mitotic doublings with 23-46% of total matings being outcrosses. These findings imply that even the mildest forms of the [PSI+], [URE3] and [PIN+] prions impart > 1% growth/survival detriment on their hosts. Our estimate of outcrossing suggests that Saccharomyces cerevisiae is far more sexual than previously thought and would therefore be more responsive to the adaptive effects of natural selection compared with a strictly asexual yeast. Further, given its large effective population size, a growth/survival detriment of > 1% for yeast prions should strongly select against prion-infected strains in wild populations of Saccharomyces cerevisiae.

Haplotype probabilities in advanced intercross populations.

Advanced intercross populations, in which multiple inbred strains are mated at random for many generations, have the advantage of greater precision of genetic mapping because of the accumulation of recombination events across the multiple generations. Related designs include heterogeneous stock and the diversity outcross population. In this article, I derive the two-locus haplotype probabilities on the autosome and X chromosome with these designs. These haplotype probabilities provide the key quantities for developing hidden Markov models for the treatment of missing genotype information. I further derive the map expansion in these populations, which is the frequency of recombination breakpoints on a random chromosome.

Gene flow from transgenic to nontransgenic soybean plants in the Cerrado region of Brazil

Evaluation of transgenic crops under field conditions is a fundamental step for the production of genetically engineered varieties. In order to determine if there is pollen dispersal from transgenic to nontransgenic soybean plants, a field release experiment was conducted in the Cerrado region of Brazil. Nontransgenic plants were cultivated in plots surrounding Roundup Ready transgenic plants carrying the cp4 epsps gene, which confers herbicide tolerance against glyphosate herbicide, and pollen dispersal was evaluated by checking for the dominant gene. The percentage of cross-pollination was calculated as a fraction of herbicide-tolerant and -nontolerant plants. The greatest amount of transgenic pollen dispersion was observed in the first row, located at one meter from the central (transgenic) plot, with a 0.52% average frequency. The frequency of pollen dispersion decreased to 0.12% in row 2, reaching 0% when the plants were up to 10 m distance from the central plot. Under these conditions pollen flow was higher for a short distance. This fact suggests that the management necessary to avoid cross-pollination from transgenic to nontransgenic plants in the seed production fields should be similar to the procedures currently utilized to produce commercial seeds.

ABSENCE OF POLLEN DISCOUNTING IN A GENOTYPE OF IPOMOEA PURPUREA EXHIBITING INCREASED SELFING.

Throughout southeastern North America, the annual morning glory Ipomoea purpurea exhibits a polymorphism at a locus that influences the intensity of floral pigmentation. Previous studies have shown that when rare, the homozygous white genotype has a greater selfing rate than the homozygous dark genotype. In the absence of pollen discounting (a reduction in transmission of pollen to other plants by genotypes that exhibit increased selfing) and inbreeding depression, this increased selfing rate should favor the white allele. Experiments reported here confirm that the white genotype has elevated selfing rates when rare but indicate pollen discounting is not associated with elevated selfing. Rather, white genotypes contribute more pollen to the outcross pollen pool. The disparity between genotypes in both selfing rates and success at pollen contribution to other plants disappears at intermediate to high frequencies of the white allele. Pollinator movements are consistent with the pattern of selfing. These results suggest that elevated selfing and enhanced success at pollen donation contribute to maintenance of the white allele in natural populations of morning glories.